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We present a semiclassical method to treat the proton breakup from a weakly bound state in 
an exotic nucleus. The Coulomb interactions between the proton, core and target are treated 
to all orders and including the full multipole expansion of the Coulomb potential. The nuclear 
proton-target interaction is also treated to all orders. The core-target interaction is included as 
an absorption. The method is semi-analytical thus allowing for a detailed understanding of the 
short range and long range effects of the interactions in the reaction dynamics. It explains also the 
origin of possible asymmetries in the core parallel momentum distributions when the full multipole 
expansion of the Coulomb potential is used. Calculations are compared to results of other, fully 
numerical, methods and to experimental data in order to establish the accuracy and reliability of 
the method. 
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I. INTRODUCTION 

The study of exotic nuclei through knockout reactions 
has provided important results along the last decades. 
Among those, one proton removal reactions have been 
proven as a useful tool for the understanding of light, 
proton rich, exotic nuclei. From the reaction point of 
view, the inclusion of a long-range part in the interaction 
between the proton and the target, i.e. the Coulomb 
field, provides a significant complication with respect to 
the neutron case as well as a challenge for reaction theory. 

This work is an extension of the method presented in 
Ref. [l[. There, a formalism to treat neutron breakup 
to all orders, including the whole multipole expansion 
of the Coulomb interaction giving the recoil of the core, 
was presented. Here that theoretical framework is com- 
pleted in order to account for proton breakup as well. 
For proton halo nuclei Coulomb breakup reactions in the 
laboratory have been used to get indirect information 
on the radiative capture, since it has been shown that 
the Coulomb breakup cross section is proportional to the 
radiative capture cross section [2j. A number of experi- 
mental papers [|[-[l3j] have reported experiments on pro- 
ton breakup from 8 B. 8 B is partner in (p,7) radiative 
capture reactions of great astrophysical interest for the 
understanding of the neutrino flux from the sun (see for 
example the discussion and references of @). 17 F is an- 
other candidate for proton halo which is still under inves- 
tigation both experimentally [l4l Il5l as well as the- 
oretically [HI, [IS US US- Furthermore 17 F plays an im- 
portant role in understanding explosive nucleo-synthesis 
in X-ray bursts and novae, as it enters in the sequence 
14 0(a,p) 17 F(p,7) 18 Ne(o;,p) 21 Na, where cross sections are 
poorly known. 

The existence of a proton halo has sometimes been 
questioned [2l| and results from different experiments 
might seem to be contradictory [22|, [23j|. Some of the 
apparent discrepancies in experimental results and their 




FIG. 1: (Colour online) Coordinate system. 



analysis have been recently discussed and in large part 
solved in Ref. [Hi. 



II. FORMALISM 

This work is the extension to the proton breakup 
problem of the formalism developed in 1] for neutron 
breakup. Similar formulae to those of the neutron 
breakup can be derived, since the Coulomb potential can 
be written as 
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where V c = Z c Z t e 2 , V v = Z v Z t e 2 and V a = (Z v + 
Z c )Z t e 2 . Pi and /3 2 are the mass ratios of proton and 
core, respectively, to that of the projectile. The coordi- 
nate system used in this paper is shown in Fig. [T] Fol- 
lowing a procedure analogous to that of [H , where x was 
defined as Xpert — \ J dte luJt V(f,t), the Coulomb phase 
for the proton is shown to be 

X p = ^-{V c e i ^ uz l v K {ub c /v)-V Q K {uR x /v) 
hv 

+V v e~^ z ^K (ujb v /v)) (2) 

with ui — {ef — £o) /h and £o is the neutron initial bound 
state energy while £/ is the final neutron-core continuum 
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energy. Since Vo = V c + V v , Eq. can be written as 

X P = x(Pi,V c ) + X (-P2,V V ) (3) 

leading to 

Xifl, V) = ^ (e^'^Koiub/v) - K (wRj_/vj) . (4) 

The Coulomb phase Eq.Q is therefore the sum of 
two terms: one of them describes the recoil of the core 
whereas the other accounts for the direct proton-target 
Coulomb interaction. Of course, in the case of the neu- 
tron the latter vanishes and the phase reduces to the one 
derived in 

The expansion of Eq. ((4]) to first order in f yields the 
dipole approximation to the phase: 

2(f3 1 V c - §2Vy) ,„ , iujz 
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which only differs from the neutron breakup case in the 
different constant factor, which is now (V c (3i — 14/32 ) in- 
stead of Vq/3! of Ref. [3. 

Our expression for the differential cross-section is 



S = 8^ I <">.\S,,U>r) + </'' 
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where |iS c t(6 c )| 2 is the core-target elastic scattering prob- 
ability discussed later. Such a form of the cross section 
has been introduced for relativistic energy electromag- 
netic excitations in Ref.(2(| where however a first or- 
der perturbation theory amplitude was considered. It 
has been derived in an eikonal approach to the halo- 
projectile scattering [U HI] and it has been widely used 
also for normal projectiles in the case of peripheral re- 
actions such as transfer and inelastic excitations [2^ |. 
Under the conditions that the core-target interaction is 
mainly absorptive it has been re-obtained by Bertulani 
in a relativistic approach to proton breakup [301 ]. This 
paper has also shown that relativistic corrections to dy- 
namical approaches such as the continuum discretised 
coupled channel method (CDCC) can increase the pro- 
ton breakup cross section of as much as 15% at the inci- 
dent energy of 250 A.MeV. Our formalism however fol- 
lows closely Ref. (2|| where relativistic kinematics is easily 
implemented and it has been used in the numerical cal- 
culations presented in this paper. As we shall see in the 
following our results agree indeed rather well with high 
energy experimental data. 

The probability amplitude in Eq.^ has been written 
as the sum of three pieces: the recoil term, 



dre- ikf VifrKe^ 



.2V C . b c 
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where, according to the discussions in p], [2^, [2?], l28l ]. 
the sudden limit has been used in order to include all 
orders in the interaction. Similarly, the second term in 
our probability amplitude is the direct proton Coulomb 
interaction. It has the same form as Eq.([7|) but for the 
substitution V c —*Vv, b c — > b v and j3i — > — 0%. 
Finally, the nuclear part is 
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A number of papers have addressed the problem of 
asymmetry in the core parallel momentum distribution 
after proton knockout 0, i, @, d M, M, HI . The fact 
that this asymmetry comes from high order terms can be 
directly extracted from our formalism. If the Coulomb 
part of the amplitude g Cou = g rec + g dlr is simply ex- 
panded to first order in x, it becomes 



g Cou ~ / dre- ik -%(?) — (V c e^ z / V K (cjb c /v) 

-V K (uR ± /v) + V v e-^ z / v K (u;b v /v)){9) 

which can be written, in terms of the one-dimensional 
Fourier transform in z— direction <bj, as 



Cou 



- V K (ujR ± /v)^(r ± ,k z ) 
+ V v K (ub v /v)<i>i(r_L, k z - p 2 u/v) 
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Thus the Coulomb breakup probability amplitude can be 
regarded as a coherent sum of three terms, each of which 
contains a shifted z— Fourier transform. The shifts are in 
opposite directions, —0\U)jv and /?2w/u, but they are not 
visible in the calculated momentum distributions shown 
in this paper as uj depends on k itself. The amount of 
deformation will also depend on the explicit form of the 
wave function 4>i{r). Moreover, the 1/v factor indicates 
that the asymmetry decreases as the beam energy in- 
creases. 

In dipole approximation, however, the amplitude is 
just 
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x I Ko(ujRx/v)-—^i(r±,k z ) 

+ K 1 (u;R ± /v)^^-4> i (r±,k z )), (11) 
R± v / 



which does not contain any asymmetry for the momen- 
tum distribution as it involves square moduli of 0i(rj_, k z ) 
separately. Hence we can confirm analytically that the 
asymmetry in Coulomb breakup parallel momentum dis- 
tributions is due to the presence of higher multipole 
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FIG. 2: Calculated inclusive momentum distribution of 7 Be 
fragments after proton-removal from 8 B against Pb at 936 
A.MeV. The calculations are the sum of contributions lead- 
ing to the ground and excited state of 7 Be weighted by the 
respective spectroscopic factors and include both the diffrac- 
tion and stripping parts of nuclear breakup, besides Coulomb 
breakup. Data are from 0. 
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FIG. 3: (Colour online) Proton momentum distribution after 
Coulomb breakup of 8 B against Pb at 936 A.MeV in both 
dipole and full-multipole approximations, for both ground and 
excited state. Spectroscopic factors are not included. 



terms, in agreement with earlier works @, 0, IE EU US H2] • 
However, the presence of the nuclear interaction intro- 
duces an interference that does depend on the sign of k z 
and thus an additional asymmetry to that due to higher 
multipole terms in the Coulomb interaction. 



III. APPLICATIONS 

This formalism has been applied to proton breakup of 
8 B against Pb and C targets at a beam energy of 936 
A.MeV, experiment reported in Q. The projectile 8 B is 
taken as a two-body object. Its ground state has com- 
ponents with a proton coupled to the ground 3/2 ~ state 
and to the excited E* = 429keV 1/2" state of 7 Be. The 
proton separation energy is 137 keV when the 7 Be is left 
in its ground state and 137 + 429 keV when it is left in 
an excited state. 

Radial wave functions have been obtained by numeri- 



cal solution of the Schrodinger equation in Woods-Saxon 
potentials with depths adjusted to reproduce the corre- 
sponding separation energies. The radius parameter of 
these Woods-Saxon potentials have been taken as 1.3 fm 
with a diffuseness of 0.6 fm. 

The nuclear breakup is described in this formalism 
through the proton-target and core-target S-matrices. 
These have been calculated by taking the optical limit in 
Glauber theory . 34j . The effective two-body interactions 
are obtained by folding effective nucleon-nucleon interac- 
tions with the densities of the nuclei involved. Pb density 
was approximated by a Woods-Saxon profile fitted to re- 
produce the experimental root mean square radius 5.5 fm 
[35| . 7 Be and 12 C densities were taken as gaussian with 
root mean square radius 2.31 and 2.35 fm respectively 
[36| . The stripping contribution to nuclear breakup has 
been calculated in the eikonal formalism (37| . 

With these ingredients and the formulae described in 
the previous section, the single particle cross-sections of 
Table [I] were calculated. Interference between Coulomb 
and elastic nuclear part are included automatically in 
the amplitude, but stripping must be added incoherently. 
When comparing to experimental data, spectroscopic fac- 
tors must be taken into account. Following the spec- 
troscopic factor for the excited state is S ex = 0.101 
in this work, while for the ground state we have used 
S gs = 0.841. The latter has been obtained by adding 
together the spectroscopic factors for two states which 
differ in core intrinsic spin Q, since our formalism can- 
not distinguish between them. These numbers provide 
an exclusive/inclusive ratio of Rtheo = 0.089 for the Pb 
target, in good agreement with the experimental result 
R 



exp 



0.085±0.021. Such agreement is also good for the 
C target, as the experimental result is R exp = 0.13±0.03, 
while we obtain Rtheo = 0.098. The inclusive momen- 
tum distribution on lead is shown in Fig. [2] with no scale 
factors. In this calculation the stripping contribution is 
assumed to have the same shape as the diffractive one. 
Clearly the agreement is not perfect, but in view of the 
very simple two-body forms of our initial and final wave 
functions which do not contain any three-body effects, it 
can be considered reasonable. Regarding absolute cross- 
sections, in the case of the Pb target, our result for the 
excited state is 57.8 mb, in agreement with the exclu- 
sive measurement of 56 ± 5 mb, while the inclusive result 
662 ± 60 mb is rightly described by our calculation of 
652.4 mb (TableHJ. For the C target we predict 9.1 (92.6) 
mb for the exclusive (inclusive) measurement, consistent 
with the experimental result 12 ± 3 (94 ± 9) mb. 

This formalism allows a direct check of the dipole ap- 
proximation for the Coulomb interaction, by compar- 
ing the results provided by the phase Eq.Q to those 
of Eq.©. Fig. [3] shows the calculated momentum dis- 
tributions of protons in the projectile reference frame af- 
ter breakup of 8 B. This does not include stripping. The 
dipole approximation, yields narrower momentum distri- 
bution and smaller single particle cross-section. High or- 
der effects, especially in the heavy target case, are not 
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TABLE I: Single particle cross-sections in mb. Beam energies 
in A.MeV. 



Projectile Target 


Energy 


&Cou+diff 




O sp 


S x cr S p 


"B 


Pb 


936 


558 


149 


707 


594.6 


8 B* 


Pb 


936 


442 


130 


572 


57.8 




C 


936 


36.9 


62.4 


99.3 


83.5 


8 B , 


C 


936 


32.1 


57.4 


89.5 


9.1 



error is of 40-60% for 17 F. Such estimates are indepen- 
dent from the initial binding energy and from the beam 
energy. Our method is an high energy approximation 
and therefore further corrections should be expected if 
the method is applied at the low energies considered in 
Ref.fH 



IV. CONCLUSIONS 



TABLE II: Barrier penetration factors wi and W2 for B and 
17 F respectively. Proton-core energies, £/, and momenta, hk, 
in MeV. ni and n 2 are dimensionless Sommerfeld parameters. 



Ef hk 


ni 


Wl 


n 2 


W2 


0.5 30.7 0.89 


0.22 


1.78 


0.013 


1.0 43.5 


0.63 


0.58 


1.26 


0.158 


1.5 53.2 


0.51 


0.81 


1.03 


0.39 


2.0 61.5 


0.44 


0.94 


0.89 


0.61 



negligible, in agreement with [38(, although in their cal- 
culation the effect appears smaller due to their choice of 
large impact parameters. 

We have also performed preliminary calculations for 
17 F. They differ from those of Ref. [3^| by about an order 
of magnitude for the g.s. breakup and by about a fac- 
tor three for the breakup from 17 F*. The reason for this 
discrepancy is in large part due to the use of a the plane 
wave as final state which neglects the Coulomb proton- 
core interaction. An estimate of the possible correction 
can be obtained by calculating a barrier penetration fac- 
tor w at the core radius R c according to Ref . flol] . 

In Table HT1 we give the values of the penetration fac- 
tors Wi and W2 for 8 B and 17 F projectiles respectively. 
We give also the values of the corresponding Sommerfeld 
parameters ni and n2, with n=Z p Z t e 2 / hv . These values 
show that for a peak in the proton-core energy spectrum 
around 1.5-2MeV there is a possible error in the cross 
sections of 10-20% in the 8 B projectile case while the 



In this paper we have introduced a semiclassical and 
analytical method to treat proton breakup to all orders 
and all multipoles in the Coulomb potential. As previ- 
ously done by some of us [l], [2?], [28| , the present method is 
based on the use of a final nucleon wave function which 
is an extention of the traditional eikonal-type. If the 
Coulomb potential is neglected, our breakup amplitude 
reduces to the well known eikonal form for the nuclear 
breakup. On the other hand if only the first order term 
of the Coulomb phase is considered, our Coulomb am- 
plitude reduces to the multipole expansion in which the 
dipole term is treated by perturbation theory, while high 
order and higher multipole terms are treated in the sud- 
den approximation according to |28| . To test the method 
we have applied it to the study of one proton breakup 
from both the ground and first excited states of 8 B. Our 
results are in very good agreement with experimental 
data at relativistic energies. Preliminary calculations for 
17 F at lower energies (65 A.MeV) have shown discrepan- 
cies of about one order of magnitude with respect to the 
results of Ref. . We believe this is mainly due to our 
choice of a plane-wave for the final state wave-function. 
However, the use of the plane- wave has allowed analyti- 
cal calculations leading to Eq. (fT0|) which is perhaps the 
most interesting result of this work. It shows explicitly 
the origin of a possible asymmetry in core parallel mo- 
mentum distributions as contained in the Fourier trans- 
form of the initial state wave function. In agreement with 
other numerical methods [19( we have demonstrated that 
the asymmetry is due to the higher order multipole terms 
in the expansion of the Coulomb potential, while it dis- 
appears in the dipole approximation. 
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